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The Morphology of the Expanding Ejecta of V2491 Cygni (2008 N.2) 
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ABSTRACT 

Determining the evolution of the ejecta morphology of novae provides valuable information 
on the shaping mechanisms in operation at early stages of the nova outburst. Understand- 
ing such mechanisms has implications for studies of shaping for example in proto-Planetary 
Nebulae. Here we perform morpho-kinematical studies of V2491 Cyg using spectral data to 
determine the likely structure of the ejecta and its relationship to the central system and shap- 
ing mechanisms. We use shape to model different morphologies and retrieve their spectra. 
These synthetic spectra are compared with observed spectra to determine the most likely mor- 
phology giving rise to them, including system inclination and expansion velocity of the nova 
ejecta. We find the best fit remnant morphology to be that of polar blobs and an equatorial ring 
with an implied inclination of 80^j 2 degrees and an maximum expansion velocity of the polar 



blobs of3100 + ^kms-' 



This inclination would 



and for the equatorial ring 2700+?q2 km s~ 
suggest that we should observe eclipses which will enable us to determine more precisely 
important parameters of the central binary. We also note that the amplitude of the outburst is 
more akin to the found in recurrent nova systems. 
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1 INTRODUCTION 

A classical nova (CN) outburst occurs in an interacting binary sys- 
tem comprising a main-sequence or evolved star (the secondary), 
which fills its Roche Lobe, and a white dwarf (WD; the primary). 
Matter transferred from the secondary to the WD builds up on the 
surface. Once extensive CNO cycle nuclear burning commences 
under these degenerate conditions, a thermonuclear runaway oc- 
curs ejecting ~ KT 5 -1(T 4 M Q of matter from the WD surface 
with velocities from hundreds t o thousands of km s (see, e.g. 
iBode & Evans 2008; Bodj20icl and references therein). A related 
sub-class of CNe are the recurrent novae (RNe), which unlike CNe, 
recur on time-scales of order tens of years. Such a short recur- 
rence time-scale is usually attributed to a high mass WD, probably 
close to the Chandrasekhar limit, togethe r with a high accretion rate 
dStarrfield et alJl985tlYaron et all2o61h . 

Nova Cygni 2008 N.2 (hereafter, V2491 Cyg) was discov- 
ered on 2008 April 10.8 UT (taken as t = ) at about 7.7 mag 
on unfiltered CCD frames dNakano et ai]|2008h and its nature con- 
firmed sp ectrometrically yielding vel oc ities of ~ 4500 km s -1 (Ha 
FWHM, lAvani & Matsumotd l2008h . iHenden & Munaril d2008h 
identified USNO-B1.0 1223-042965 as the likely progenitor sys- 



te m (within 0.9 arc seconds of V2 491 Cyg). An archival search 
by ljurdana-Sepic & Munarf d2008l) of the Asiago Schmidt tele- 
scope images of the progenito r collected over the p eriod 1970-1986 
showed no previous outburst. Munari et al. (2010b) have provided 
a detailed study of the photometric and spectrometric evolution of 
V2491 Cyg, including a photo-ionization analysis of the ejecta and 
their chemical abundances. They derived the reddening for V2491 
Cyg of E(B -V) = 0.23±0.01, a distance of 14 kpc, and V = 17.9, 
R c = 17.5, I c = 17.1 for the progenitor in quiescence. 

In Fig. [T] we show the light curve for the first 50 days af- 
ter outburst from the AAVSQj. A rapid decline is observed for 
the first 10 days after outburst and then a sudden rebrightening 
peaking around day 15. It is not fully u nderstood why this re- 
brightening should happen. For example, Hachi su & Katol J2009h 
have modelled the rebrightening by introducing magnetic activ- 
ity as an additional energy source to nuclear burning. Although, 
a magn etic WD has a s pin closely related to the binary orbital 
period dKing et aT] 1990l). taken p reviously by several authors as 
0.0958 days dBaklanov et alj|2008l) , but there is no observed evi- 
denc e of a short period i city in the X-ra ys related to the orbital pe- 
riod ( Ibarra et al . 2009; Pag e et alj|20ld) and generally polars are 
weak X-ray sources dicing & Watsonl ll987). It is also noteworthy 
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Figure 1. Optical light curve for V2491 Cyg from AAVSO data. The verti- 
cal lines are days where data were taken with the Meaburn Spectrograph on 
the 2-meter robotic Liverpool Telescope. 



that A. Baklanov (private communication) is not confident that the 
above period is rel ated to the orbital pe riod of the system, a matter 
that is addressed in lDarnlev et al.l d20ld) . 

V2491 Cyg is of particu lar interest because it was detected as 
an X-ray sou rce pre-outburst dlbarra & Kuulkersl200 8 ; Ibarr a et al.l 
I2008L l2009h . This is only the second C N to be detected in X - 
rays pre-outburst (after V2487 Ophiuchi. iHernanz & Salall2002h . 
V2487 Oph was originally suggested to be a RN because of its 
rapid decline in the optical d omain and the presence of a plateau 
phase during decline (see e.g jHachisu et alj200 2; Hernanz &Salal 
2002). V2487 Oph's true recurrent nature came to light only after a 
search of the Harvard College Observatory archiva l photographic 
collec tion revealed an outburst on 1900 June 20 dPagnotta et al.l 
2009). They used two methods to determine a recurrence timescale 
of order ~ 20 years, a Monte Carlo simulation to calculate the 
probability that a given recurrence timescale would produce ex- 
actly two discovered eruptions and by directly estimating the most 
likely number of eruptions. Furthermore, V2491 Cyg has very sim- 
ilar proprieties to the suggested extragalactic RN M31N 2007- 12b, 
which was shown to have similar magnitud e and colour to the 
RN RS Oph at qu iescence dBode et al]|2009h . On the other hand, 
Ibarra et al. I d2009» derived 0.2-10 keV X-ray luminosities rang- 
ing from 10 34 -4xl0 35 erg s _1 implying inter-outburst mass ac- 
cretion rates in the range of 10~ 9 to 10 -8 M Q yr~' for a 1 M Q 



WD. Such low values for the accretion rate would imply a recur- 
rence timescale of > 100 yrs dYaron et alj|2005l ; llbarra et alJI 20091 ; 
iPage et alj[20lch . However, a higher WD mass will reduce the re- 
currence time. 

Post-outburst Swift observations of V2491 Cyg showed a 
clearly detected X-ray source, although the count rate was almost 
an order of magni tude fainter than during pre-outburst observations 
dPage et al.l 2010). The post-outburst observations showed three 
distinct phases: first a hard X-ray emission phase, then a super- 
soft source (SSS) phase and finally the decline of the SSS with the 



hard component again becoming prominent, since it fades much 
less rapidly than the SSS emission dPageetal.ll2010l) . 

For decades it has been known that spectral line profiles con- 
tain a great deal of information abou t the expanding remnant of a 
nova outburst (e.g., lHutchings|[l972T) . Here we use spectrometric 
data to constrain the structure and velocity field for the expanding 
remnant of V2491 Cyg. In Section[2] we describe our observations 
and data reduction methodology. The results are presented in Sec- 
tion [3] and then modelled in Section [4] In Section [5] we present the 
results of the modelling and in Section [6] we discuss the implica- 
tions of the results. 



2 OBSERVATIONS AND DATA REDUCTION 

Optical spectra of V2491 Cyg were obtained with the Meaburn 
Spectrograph (a thermocooled, Apogee AP7p 512x512 pixel 
thinned Tektronix array prototype fibre-fed robotic spectrograph, 
now superseded by FROD OSpec) on t he 2-metre robotic Liver- 
pool Telescope (LT; ISteele et al]|2004h sited at the Observatorio 
del Roque de Los Muchachos on the Canary Island of La Palma, 
Spain. The spectrograph was fed using a close-packed fibre bun- 
dle input array, consisting of 49 x 0.8" diameter fibres which were 
reformatted as a slit with the fibres in random order. 

The observations were carried out over several epochs from t 
= 5 days (2008 April 16) to 31 days (2008 May 12) post-outburst. 
The log of spectrometric observations is shown in Table [T] The 
spectrograph provided three overlapping fixed grating positions, 
giving complete wavelength coverage from 3900-8000A. How- 
ever, we only use a single grating position (5350-6888A) with dis- 
persion of 3.00 A/pix and spectral resolving power ~ 350. Data 
reduction was performed through a pipeline that initially performs 
bias, dark frame and flat field subtraction. With a comparison arc 
spectrum, a third order polynomial 2D wavelength/distortion map 
was created using the STARLINK Figaro ARC/IARC/ISCRUNCH 
routines and applied to the object image. Sky subtraction based 
on identification of non-object areas of the 2D spectrum was then 
carried out, if significant sky lines were visible, using the Figaro 
routine POLYSKY. A ID spectrum was then automatically ex- 
tracted using the PROFILE and OPTEXTRACT routines. Rela- 
tive flux corrections were accomplished by using observations of 
a number of early type Be stars (showing few absorption lines) ob- 
tained for a different programme. Absorption lines were patched 
over and t he star's temperature w as chosen following their spec- 
tral types ( Boeh m-Vitensd 1 198 ll) . Using the STARLINK Figaro 
IRFLUX routine we corrected the spectra for instrumental effi- 
ciency/atmospheric absorption. In Fig. [2] we show the computed 
flux-corrected spectra from t = 5 to 31 days after outburst. 

We also obtained spectrometric observations at t = 2, 25, 33, 
108 and 477 days after outburst with: (/) the 0.6m telescope of 
the Schiaparelli observatory in Varese equipped with a multi mode 
spectrograph and various reflection gratings, (ii) the 3.5m TNG 
in La Palma and the high resolution spectrograph SARG (AA/A 
= 75000), and (Hi) the spectrograph/imager AFOSC on the 1.82m 
telescope of the Padova Astronomical Observatory in Asiago with 
wavelength coverage of 6400-7050A. This is part of the obser- 
vationa l monitoring of V2491 Cyg is described by iMunari et al.l 
d2010bl) to which the reader is addressed for full details about data 
acquisition and reduction. 
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Figure 2. Meaburn spectrometric observations of V2491 Cyg during days 5-31 (2008 April 06 to 2008 May 12). The Meaburn spectra have been corrected 
for instrumental efficiency/atmospheric transmission using Be Star comparison spectra (see text for details). On the left-hand-side we have truncated the Ho- 
emission to show more clearly the lower intensity lines. The spectrum on day 31 after outburst was too noisy at this low level, so we only plot the Ha line. 
On the right-hand-side the Ha line profiles are shown. The time since outburst (in days) is shown to the left of each spectrum. Cosmic rays have also been 
removed. The higher resolution (non-Meaburn spectrograph) Ho- profiles in Table[T]are shown in Fig. [4] 



Table 1. Log of optical spectral observations of V2491 Cyg. 



Date of 
observation 


days after 
outburst 


Exposure 
time (sec) 


Instrument 


2008 April 13 


2 


900.0 


Varese-MMS 


2008 April 16 


5 


200.0 


Meaburn-LT 


2008 April 19 


8 


200.0 


Meaburn-LT 


2008 April 20 


9 


200.0 


Meaburn-LT 


2008 April 21 


10 


200.0 


Meaburn-LT 


2008 April 22 


11 


200.0 


Meaburn-LT 


2008 April 23 


12 


200.0 


Meaburn-LT 


2008 April 24 


13 


200.0 


Meaburn-LT 


2008 April 25 


14 


300.0 


SARG-TNG 


2008 May 03 


22 


600.0 


Meaburn-LT 


2008 May 06 


25 


600.0 


Meaburn-LT 


2008 May 12 


31 


600.0 


Meaburn-LT 


2008 May 14 


33 


1800.0 


Varese-MMS 


2008 July 27 


108 


900.0 


Asiago-AFOSC 



3 OBSERVATIONAL RESULTS 

We analysed the spectral information first to get an initial estimate 
of the component parameters before detailed modelling of V2491 
Cyg began. The spectra show broad Ho-, nitrogen (Nn 5679A) 
and sodium ( Nai D ) lines (Fig. |2j. The spectral monitoring by 
Munari et aDHpiOb) shows that at day +33 after outburst the au- 
roral [Nn] 5755 line was still pretty weak, and thus the nebular 
[Nn] 6458, 6584 lines should not yet be perturbing the Ha pro- 
file. Also evident is a narrowing of the Ho- line profile from day 
22 after outburst, which is after the optical rebrightening. We can- 
not say why this rebrighte ning should happen but it is noteworthy 
that iMunari et al.l d2010bl) also observed a decrease in the equiv- 
alent widths of the H/J line around this time. The Nn and Nai 
lines seen are similar to th ose seen in the early spectrum of V394 
CrA fwilliams e t al.|[l99lh . However, these lines disappear from 
the spectrum of V394 CrA 20 days after its 1987 outburst. In the 
case of V249 1 Cyg both are clearly present later after outburst. The 
N n line, as in V394 CrA, is most likely due to continuum fluores- 
cence from the post-outburst WD radiation. 

The line of interest for our detailed analyses is Ha because of 
its strength compared to those of other lines visible in the spectrum. 
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Table 2. Derived radial velocity for each of the 5 components of the Ha line in V2491 Cyg (see Fig.[3]and text for details). 





days after 




Component Radial Velocity (km s ) 




outburst 


1 


2 


3 


4 


5 


2008 April 16 


5 


-1875±32 


-530±34 


117+3529 


615±64 


2028±68 


2008 April 19 


8 


-1860+17 


-510±22 


112±961 


671+61 


2085±19 


2008 April 20 


9 


-1899±50 


-537±38 


98±849 


700±28 


2()88±3 


2008 April 21 


10 


-1828±65 


-482±48 


99±2529 


653±8 


1994±53 


2008 April 22 


11 


-1909±76 


-531 ±55 


92±1209 


687±59 


2043±24 


2008 April 23 


12 


-2025±163 


-642±93 


120±1551 


680±87 


2072±120 


2008 April 24 


13 


-2054±67 


-720±72 


90±1341 


558±29 


2009±65 


2008 May 03 


22 


-1634±47 


-610±30 


184±29 


673±30 


2088±26 


2008 May 06 


25 


-1739±218 


-716±20 


96±33 


687±58 


2015±44 


2008 May 12 


31 


-1865±92 


-822+111 


334+91 


368±368 


2094±5 



Further more, later in outburs t other emission lines are severely 
blended dMunari et ai]|2010bl) . Using the low resolution Meaburn 
spectra we initially fit five Gaussian components to the spectra (see 
Fig. O. This allows us to derive the individual components' ra- 
dial velocities and their respective FWHM (see Tables [2] and [3] re- 
spectively). There is evidence for velocity symmetry between the 
faster moving components, 1 and 5. Components 2, 3 and 4 are 
much harder to constrain, especially component 3 which becomes 
stronger after rebrightening. We find the displacement of the line 
profile from the rest wavelength of the Ha, using components 1 
and 5, to be 92±36 km s~' . 



4 MODELLING 



Using shape (Version 3.56, Isteffen & Loped 120061 : ISteffen et all 
201013 we performed morpho-kinematical studies of V2491 Cyg. 
Class ical nova ejecta have been modelled with various structures 
(e.g., iHutchingslI 19721 : [So3l 19831 : iGill & 0'BrierJ[l99^) . Work by 
ISlavin et alj h995l) . later updated bv lBodd d2002h . suggested a re- 



lationship between the speed class (in terms of the time for the 
wind flux to decline 3 magnitudes from peak - fa) and the major to 
minor axis ratio of the expanding nova shell, where the faster the 
nova the less the degree of shaping. This relationship works well 
with CNe (short orbital period systems) and possibly RNe of short 
orbital peri od. However, RS Oph, a system with orbital period of 
455 days dKenvorJ[l994l : iFekel et ai]|2000l : iBrandi et alj|2009l) , at 
155 days afte r outburst showed a deprojected major to minor axial 
ratio of 3.85 teibeiro et all 2009). which considering its fast opti- 
cal dechne_wouldjiot agree with the relationship originally found 
bv lsiavin et al.l dl995h . Recent res ults for V2672 Oph also appear 
to disagree with this relationship ( Mu nari et al. I l2010ah . However, 
we expect that the mechanism for formation of CN and RN shells, 
especially those RNe with long orbital periods, to be different. For 
example, systems such as RS Oph have remnants grossly affected 
by the interaction of ejecta with the pre-existing red-giant wind, 
which is unique to this sub-group of RNe. 

Several mechanisms have been proposed to be involved in 
the shaping of CN ejecta. These can be identified as that associ- 
ated with the common-envelope phase during outburst, the pres- 
ence of a magnetised WD, the spin of the WD and an asym met- 
ric thermonuclear runaway (see e.g., lO'Brien & Boddl2008l) . The 
most widely accepted model involves the common-envelope phase 



Available from http://bufadora.astrosen.unam.mx/shape/ 
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Figure 3. Gaussian components best fit to the observed Hey line profile for 
V2491 Cyg on day 25 after outburst. The observed spectrum is shown in 
black and the combined Gaussian components in dashed green. The colour 
scheme for the individual Gaussian components are just to separate the dif- 
ferent components. The component numbers, 1 through 5, are in increasing 
wavelength of the peak. Also shown are the residuals to the overall fit. The 
results for the derived radial velocities of each peak and FWHM of each 
component are presented in Tables|2]and|3] respectively. (A colour version 
of this figure is available in the online journal.) 



where the ejecta from the WD envelope engulf the secondary star 
within a matter of minutes following outburst. Th e secondary trans- 
fers energy and angular momentum to the ejecta dLivio et aljl990l : 
iLlovd et al]| 19971 : IPorter et"al]|l998l) . WD rotation, when incorpo- 
rated into calculations of the common-envelope phase, produces 
the observed prolate remnants jPorter et aljl 998). 



4.1 Model Assumptions 

It is desirable when dealing with modelling to reduce computing 
time and yet have a physically sound model. The modelling of per- 
mitted line profiles so early after outburst is complicated by the fact 
that they are most likely affected by an optically thick medium. The 
full treatment of radiative transfer requires full knowledge of sev- 
eral parameters of the system and further assumptions are required. 
For example, the temperature of the central system in outburst, 
which can be derived from multi-frequency observations, and the 
density of each component, which is much harder to derive since 
we must assume an ejected mass and a geometry of the system. 

We considered the effects of radiative transfer on the line pro- 
file and found that the the lower the inclination the more the ab- 
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Table 3. Derived FWHM velocity for each of the 5 components of the Ho- line in V2491 Cyg (see Fig.[3]and text for details). 



days after Component FWHM Velocity (km s ) 

outburst 12 3 4 



2008 April 16 


5 


1729±47 


1357±24 


1870±3487 


1379±130 


1166+91 


2008 April 19 


8 


1664±48 


1412+14 


1807±1780 


1429±109 


1115±41 


2008 April 20 


9 


1649±40 


1399+101 


1998+1222 


1488±57 


1050±27 


2008 April 21 


10 


1658±101 


1239+53 


1995±2238 


1261+137 


1 194±74 


2008 April 22 


11 


1660±89 


1283+59 


2694±2978 


1305+137 


1084±53 


2008 April 23 


12 


1639±180 


1384±210 


1617±776 


1476±395 


1030±167 


2008 April 24 


13 


1421±110 


1335+83 


1765±2077 


1364+169 


1127±92 


2008 May 03 


22 


1340±60 


910±27 


931 ±86 


1502±105 


806±40 


2008 May 06 


25 


1232±255 


996±236 


1081±678 


1477±382 


875±63 


2008 May 12 


31 


1189±103 


1096±76 


1725±298 


1832±324 


795±26 




6450 6500 6550 6600 6650 6700 

Wavelength (A) 



Figure 4. Spectral evolution of the high resolution Ha line profile from 
several telescopes as shown in Table[T]and in the text. The spectra has been 
offset and the numbers on the left are days after outburst. 

sorption affects the system but the overall line shape does not sig- 
nificantly change from the optically thin case. Furthermore, the ab- 
sorption affects more the lower velocity systems. So when dealing 
with the results caution is necessary especially with the densities 
derived from the models because applying a density in the opti- 
cally thin case will actually be only a lower limit to the true density 
of a system with greater opacity. 

We also made the assumption that at ear ly times the structures 
are still volume filled, della Valle et al. have shown the Ha 

spectral evolution of V382 Vel from early outburst to the nebular 
stages (5 to 498 days after outburst); at early stages the line profile 
showed an asymmetric saddle-shape with the blue component more 
prominent than the red. The blue component diminishes in strength 
and the whole line profile develops to become a flat-topped shape. 
They associate this with the termination of the post-outburst wind 
phase and the complete ejection of the envelope. Our higher resolu- 
tion spectra, Fig. [4] also show a similar development, and a similar 
interpretation can be made here. 

4.2 Model Parameters 

Given V2491 Cyg's fast decline from maximum we assume an 
axial ratio o f 1.15 for o u r cho sen structures following the rela- 
tionships in ISlavin et alj |l995). We also consider the possibil- 
ity of a structure similar to RS Oph, first coined "peanut-shaped" 



a) b) 




Figure 5. Models as visualised in shape, a) polar blobs with an equato- 
rial ring (Model A), b) a dumbbell structure with an hour-glass overdensity 
(Model B), c) prolate shell with an equatorial ring (Model C) and d) prolate 
shell with tropical rings (Model D). 

(Bod e et al 1l2007h and later modelled as a dumbbell structure with 
an hour-glass-shaped central overdensity jRibeiro et alj2009f) . This 
model is included in our studies because, as mentioned before, this 
system has been suggested as a R N and a similar structur e has 
been proposed for the RN U Sco jDrake & Orlandoll201(]h . The 
candidate structures we chose are (i ) polar blobs with an equato - 
rial ring (Model A, e.g., HR Del; iHutchingsl fl972l ; ISoldll983h : 
(ii) a dumbbell with an inner hour-glass overdensity with an ax- 
ial ratio the same as derived above for RS Oph (Model B); (iii) 
prolate shell with an equatorial ring (Model C); a nd (iv) a prolate 
shell with tropical rings (Model D, e.g., V705 Cas; iGill & O'Brienl 
1999). As noted above, Model B differs from that generated by 
the common-envelope phase described above because the observed 
structure could arise due to interaction with a pre-e xisting red- 
giant wind and/or collimation by the accretion dis c (e.g. jBode et al.l 
2007; Sokoloski et al. 20081; iRibeiro et alj|2009h . Our models, as 
visualised in shape, are shown in Fig. [5] 

The lack of a resolved image of the remnant means any model 
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is not as well constrained as when the structure is spatially resolved. 
However, for the models assumed above we constrain a system in- 
clination (where an inclination i = 90 degrees corresponds to the 
orbital plane being edge-on, and i = degrees being face-on) and 
an expansion assumed to take place in a Hubble flow. 

We explore the parameter space from 0-90 degrees and max- 
imum expansion velocity (Vexp) ranging from 100-8000 km s 
(in steps of 1 degree and 100 km s -1 , respectively). We exclude 
Vexp > 8000 km s _1 as this is much higher than those normally 
associated with novae and indeed is a t the lower end of veloci- 
ties found in Type la Supernovae (e.g., lLividl2000h . We note here 
however, that in the helium nova V445 Pup high speed knot s were 
observed moving at 8450±570 km s" 1 dWoudt et al]|2009l) . The 
models were run several times to produce a well sampled model 
spectrum. The model spectra are compared to the observed spec- 
tra and flux matc hed via X 1 minimisation, using techniques from 
IPress etaljjl992h . Furthermore, an estimate of the density ratio of 
each component was derived from their total fluxes using the task 
specfit in These were used as inputs for the relative struc- 

ture densities. 



5 MODELLING RESULTS 

5.1 First epoch observations at r = 25 days after outburst 

Shown in Fig. [6] are the results of our model fits, and their respec- 
tive lcr errors, for all the model runs. Panel a) shows the results for 
Model A. These suggest that the best fit for this model is an incli- 
nation of 80+12 degrees with Vexp for the polar blobs of ~ 31 00+™ 
km s" 1 and for the equatorial ring of ~ 2700+™ km s -1 . The in- 
dividual component fluxes are about the same. Panel b) shows the 
results for Model B. These suggest that the best fit to the model is 
for a system with inclination 63+J 3 degrees and V e xp of 3500++3g[j 
km s _1 . The flux is dominated by the dumbbell component. Panel 
c) shows the results for Model C. These suggest that the best fit to 
the model is for a system with inclination 50+^ degrees and V e xp 
of 2800+2™ km s -1 . The flux is dominated by the prolate structure. 
Finally, panel d) shows the results for Model D. These suggest that 
the best fit to the model is for a system with inclination 54+2* de- 
grees and Vexp of 2800+2™ ' an 8 • ^ ne " ux K a 8 am dominated by 
the prolate structure. 

As previously mentioned these results are not as well con- 
strained as they would be if we could have the spectra combined 
with imaging data. We also performed modelling for a simple 
sphere. The results suggest that the best fit maximum expansion ve- 
locity Vexp ~ 2700 km s -1 . The shape of the synthetic spectra was 
that of an inverted U (similar to the prolate components in Models 
C and D of Fig.©. 

To determine the best fitting morphology we examined the 
derived reduced X 1 probabilities and found that Models A and B 
showed comparable probabilities. The other structures gave negli- 
gible probabilities. 

5.2 Second epoch observations at r = 108 days after outburst 

Another test we perform on our models to see which best repro- 
duces the data, considering only Models A and B, is to evolve 

3 IRAF is distributed by the National Optical Astronomical Observatories, 
which is operated by the Associated Universities for Research in Astron- 
omy, Inc., under contract to the National Science Foundation. 



the models to a later epoch, keeping the parameters for inclina- 
tion and Vexp as derived from the first epoch to the second epoch. 
In th e case of Model A, it was allowed to evolve linearly. How- 
ever, iRibeiro et alj ( 120091) suggested that their central component 
for RS Oph showed deceleration and we therefore kept this central 
component the same size from one epoch to the next in Model B. 

What is evident is that if the models are allowed to evolve lin- 
early, with the derived parameters the same, then the line profile 
should look identical from one epoch to the next in the absence 
of other physical changes. For the line profile to change, but with 
the velocity field and inclination constant, there must be a transi- 
tion involving other factors. Just changing the densities does not 
suffice without the addition of new components. This transition 
may therefore be associated to the termination of the post-outburst 
wind phase and the complete ejection of the envelope, which is 
propos ed to happen within the first couple of weeks (as in V382 
Vel; Idella Valle et al.l l2002). Therefore, the later epoch modelling 
does not have a filled distribution but now has a depth of mate- 
rial (shell) determined by the first epoch model. We therefore apply 
a volume distribution within this shell to both models and com- 
pare with the spectrum at d ay 108 after outburst (Fig.[4]l. Further- 
more, |Mjmarj_et_al. ( 2010b) from photo-ionisation models showed 
the nova to be in the nebular stage at day 108 after outburst and the 
[Nil] 6548 and 65 84 A emission lines may contribute significantly 
to the Ha profile. Certainly [O m] was about ten tim es stronger than 
Ha at this time (see Table 4 in Munari et al. 201 Obi). 



5.2.1 Model B - Dumbbell with an hour-glass overdensity 

We do not show the Model B synthetic spectrum as it did not repro- 
duce the observed second epoch spectrum at all well. The synthetic 
spectrum was that of an upside down "T" shape. We therefore, be- 
lieve that Model A is the correct model for V2491 Cyg. 



5.2.2 Model A - Polar blobs and equatorial ring 

In light of the presence of [N n] we should account for this when 
modelling the spectra on day 108. We should expect [Nn] to appear 
as the expanding remant becomes less dense overall and the nebular 
spectrum emerges. In the first instance, to fit the observed spectrum 
the velocities were kept constant from the first epoch for the Ha 
and [Nil] lines and only the relative flux was allowed to change. 
We assume here that the model derived for Ha is the same for [N n] . 
This did not reproduce the observed spectrum without introducing 
significant additional absorption features. 

We therefore allowed the velocity of the different components 
to change keeping the flux ratio of the [Nn] lines always at 1:3 and 
the relative flux of the Ha to the [Nn] lines as a free parameter. 
The result of this modelling is shown in Fig. [7] This provides an 
excellent fit to the spectrum at day 108. The derived velocities for 
the polar blobs were ~ 3500 km s" 1 (Ha) and ~ 3600 km s _1 ([Nn]) 
and for the equatorial ring ~ 2700 km s _I (Ha) and ~ 2600 km s _I 
([Nil]). These results are within 2<r of those derived for the first 
epoch. 

It should be noted that there is a change in the contributions 
to the flux of the different components of Ha. Where in the first 
epoch the two components were of comparable flux, in the second 
epoch the ring is the major contributor to the flux with relative flux 
contribution of the polar blobs to the ring of 1:2.5 for Ha. This 
change in relative flux can be easily understood as arising from 
the fact that material moving at higher velocities will fade faster 
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Figure 6. Results of the best fit model spectrum to the observed spectrum on day 25 after outburst. The images are for a) polar blobs with an equatorial 
ring (Model A), b) dumbbell structure with an hour-glass overdensity (Model B), c) prolate shell with an equatorial ring (Model C) and d) prolate shell with 
tropical rings (Model D). Displays, top - the most likely result for the structure, where the colour gradient represents the probability that the observed X 1 
value is correct. Bottom - the observed (solid black) and modelled spectrum (short-dash black). Also shown are the individual component contributions; rings 
(dot-dash green), prolate (long-dash blue) and polar blobs and dumbbell (short-dash-long-dash red). (A colour version of this figure is available in the online 
journal.) 
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Figure 7. Fit to the observed line profile (grey) at day 108 after outburst. The model spectrum (black) is the sum of the individual components, Ha (green), 
[Nn] 6584A (red) and [Nn] 6548A (blue). The ring (dashed line) and polar blobs (solid line). The ratio of the [Nil] lines was kept at 1:3 and the relative flux 
of Ha to the [N n] lines was kept as a free parameter. (A colour version of this figure is available in the online journal.) 



in Ho- emission. Therefore, the change shown in Fig. [7] where the 
equatorial ring becomes the dominant component stems from the 
fact that it presents an expansion velocity lower than that of the 
polar blobs. 

In the case of the [Nn] line the polar blobs to ring relative 
flux contribution is 2.2:1. Their expansion velocities are ~ 3600 
km s -1 and ~ 2600 k m s -1 , respectively. We should also note that 
iMunari et al] feOlObL their Figure 10) show a line profile at day 
477 which is very different from the day 108 line profile. It does 
not have a flat top profile and has significant asymmetry. Our model 
can also account for this if the densities of the different components 
continue to evolve. As the nebula becomes less dense we expect 
that the relative contribution from [N n] should increase. 

There is one more test we can perform to further constrain our 
model. The derived inclination of 80 degrees is very high and we 
therefore expect eclipses whose existence we are now investigating. 
Such eclipses should be easily observed in this system. For exam- 
ple, the CN DQ H er (orbital period ~ 4 h) shows deep ~ 0.9 mag 
eclipses in UBVR dZhang et al.ll 19951) as does the RN U Sco with 
eclipses of depth ~ 1.5 mag in B dSchaefer & Ringwal j[l995l) . Of 
course, if V249 1 Cyg turns out to be eclipsing this will also enable 
us to determine more precisely important parameters of the central 
binary. 



6 DISCUSSION AND CONCLUSIONS 

Of the morphologies discussed above early in outburst, Models A 
and B, are the best fitting and give results that are very similar to 
each other early in outburst and this is reflected on the fact that 
the inclinations of the two morphologies overlap, within the errors. 
Furthermore, the main real difference is in the size of the central 
region. Model A is representative of a CN system, while the Model 
B has a morphology that has been observed in the RN RS Oph. 

To model the second epoch observations the model was al- 
lowed to evolve into a system that was not volume filling but a 
shell. To replicate the observed spectrum with Model A at this 
time we needed to introduce significant emission from the [Nn] 
lines and the components densities required to be changed with the 
equatorial ring becomimg the dominant contributor to the flux. The 
change in derived densities between the components should also 
be expected since the different components are moving at differ- 
ent speeds. Emission from the polar blobs moving at higher speed 
will be observed to fade faster than the equatorial ring. Model B at 
the second epoch did not reproduce the line profile. Therefore, from 
the second epoch observations we suggest that Model A is the most 
likely with an inclination of 80^, degrees and a maximum expan- 
sion velocity of 3100^™ km s , for the polar blobs. This velocity 



The Morphology of the Expanding Ejecta ofV2491 Cygni (2008 N.2) 9 



is in line with t hose found in hi gh-ionisation absorption lines in 
X-ray spectra bv lNess et alj ( 120101) . 

The high inclination found for Model A is however difficult 
to reconcile with the ou tburst amplitu de vs. decline rate relation- 
ship for classical novae dWarnerll 19871) . The observed amplitude of 
10 mags is much lower than that predicted for the derived incli- 
nation of 80 degrees (15.5 mags) while a pole-on system would 
still have a predicted amplitude of 12.5 mags. However, we note 
that RNe p resent a similar o utburst amplitude at a similar inclina- 
tion (e.g. MunarietalJl999h . This is in line with the suggested RN 
nature of this object from other authors (e.g., Page et al. 2010]). Fur- 
ther exploration of th e nature of the progenitor will be presented in 
|Parnlevetal] ( l2010h . 
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